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Initial states of high energy heavy ion collisions are studied using a dipole model through the 
DIPSY event generator that dynamically includes saturation together with the fluctuations and 
correlations of the BFKL cascade. The eccentricities £1,2,3,4 are calculated at RHIC and LHC. 
Predictions are made for correlations and fluctuations in rapidity of the eccentricities, and conven- 
tional theoretical approximations are tested. A large set of initial state Au-Au, Cu-Au and Pb-Pb 
collisions have been generated and are published online. 



I. INTRODUCTION 

The study of azimuthal correlations and their origin 
plays a big role in high energy heavy ion collisions. Often 
the azimuthal correlations are explained by an anisotropy 
in the transverse plane of the state just after collision, 
that propagates to the final state through collective ef- 
fects ft 

Traditionally these studies have been made with av- 
erage transverse shapes as function of impact parameter 
b, but lately more and more studies have incorporated 
event- by-event fluctuations of the interaction shape 
The first sign was that fluctuations are necessary to ex- 
plain the large elliptic flow in central Cu-Cu collisions 
seen at PHOBOS [J], but recently also fluctuation driven 
observables such as triangular flow [5l-fl0l| and directed 
flow [HI EH have been studied. Some of the approaches 
to model these fluctuations are: 

The colour glass condensate model fl3| . for example 
implemented in the KLN model calculates quanti- 
ties from an average transverse gluon density, evolved 
through rapidity. This approach does not include the 
event- by- event and region- by- region fluctuations and cor- 
relations in BFKL. 

NeXus [l5[ is an initial state event generator based on 
a Gribov-Regge model of hadronic collisions, and is now 
combined with the hydrodynamical model SPheRIO [l6[ 
to make the full final state event generator NeXSPhe- 
RIO [13 ■ UrQMD Ultra-relativistic Quantum Molecu- 
lar Dynamics [a, LL8] uses flux tube excitation and frag- 
mentation to generate initial states. Another model is 
AMPT fig. |20I |. a multi-phase transport model, using 
HIJING [21j for initial states, where nucleon-nucleon col- 
lisions are found from an eikonal formalism. 

NeXus, UrQMD and HIJING all use binary collision 
between the nucleons, and thus do not take saturation in 
the cascade into account properly. 

In this letter the dipole model introduced in [23 - |25| . 
with the Monte Carlo implementation DIPSY, will be used 
to describe the states at time t = 0. This dipole model is 
based on BFKL in impact parameter space, and includes 
all the fluctuations and correlations from a partonic cas- 
cade while still taking saturation into account in both 
evolution and interaction. 

The model is briefly introduced in section|TT]and results 



of common quantities such as eccentricities £1,2,3,4 follow 
in section IIIII In section IIV1 the effects of saturation on 
the BFKL correlations and fluctuations are studied, and 
section fVl ends the paper with concluding remarks. 

II. THE MODEL 

The model and its Monte Carlo implementation DIPSY 
are described in detail in [22|]. Below follows a short 
summary of the key components, and why it is well suited 
for the analysis of azimuthal eccentricities in high energy 
heavy ion collisions. 

DIPSY is an event generator using colour dipoles in 
transverse space. The two starting dipole states, repre- 
senting the colliding particles, evolve in rapidity towards 
each other to meet at the interaction rapidity yo. The 
interacting dipoles in the virtual cascade are traced back 
towards their valence parents, separating the real gluons, 
referred to as backbone gluons, from virtual fluctuations. 

The evolution and interaction amplitudes are based 
on leading logarithm BFKL, and is corrected for several 
nonlcading effects: 

• Running a s . The scale is set by the dipole sizes. 

• Energy conservation. The full 4-momentum and 
recoils of each parton are tracked, and p+ and p- 
ordering is required, simulating non-singular terms 
in the splitting function and "energy scale terms" . 

• Confinement is included through a gluon mass, sup- 
pressing large dipole emissions. 

• Saturation. Unitarisation provide saturation in the 
interaction, and the "dipole swing" allows for merg- 
ing of gluon chains in the cascade. 

The running a s and the lightcone momentum ordering 
makes the model take the essential parts of non-leading 
logarithmic effects into account, while still including the 
fluctuations from a BFKL description. Multiple interac- 
tions together with the dipole swing gives a description 
of merging and splitting gluon chains at any rapidity, 
independently of where the interaction rapidity j/o is set. 

A proton in DIPSY starts out as a triangle of dipoles, 
that then evolves in rapidity to form more complicated 
states before collision. A heavy ion is here described 
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FIG. 1: (Colour online) The eccentricities at Au-Au collisions 
for ^/sjvjv = 200 GeV as function of impact parameter b. The 



error bars show ov 



VW) - <^> 2 



FIG. 2: (Colour online) The eccentricities at Pb-Pb collisions 
for y^ijvjv = 2.76 TeV as function of impact parameter b. The 
error bars show a En = yj (e|) — {£n) 2 - 



by a Wood-Saxon distribution of nucleons with a hard 
core [26|, [27j , where each nucleon is described by a tri- 
angle just like a proton. This set of triangles in impact 
parameter space is then evolved as one state giving a web 
of gluon ladders and loops connecting the nucleons. 

The backbone gluons describe the t = state, which 
will be used to study properties such as eccentricities 
in this letter. The gluonic states can also be used 
as initial condition in a final state evolution model. 
For this purpose, a large set of collisions of several 
reactions has been generated and published online at 
http://home.thep.lu.se/ christof/DIPSYEvents/ 



III. RESULTS 

A. e n at RHIC and LHC 

The calculation of e n in this letter will follow the 
method in [f| with the formula 

yj (r 2 cos(ri(/))) 2 + (r 2 sin(n(j))) 2 



Here the averages are over all real gluons in a rapidity 
slice 77 G [—1,1], and r and cj) for the gluons are deter- 
mined with respect to the center of gravity. 

The eccentricities will be shown as function of the num- 
ber of participants. The spectators can be identified in 
DIPSY as the nucleons that have no interacting emissions. 

Also the formula 

yj (r n cos(n0)) 2 + (r n sin(n0)) 2 

is in use [6, 28] for the eccentricities. The difference be- 
tween the two definition is essentially only a constant 



factor (n + 2)/4 as argued in [28|, and this ratio is con- 
firmed by DIPSY data. 

£1 however, is calculated according to [ill, EH from 

yj (r 3 cos(0)) 2 + (r 3 sin(0)) 2 



For a closer comparison to observables, yj (e 2 ) will be 
shown rather than (e„) [28|. This adds about 15%, 7%, 
13% and 10% to £1,2,3,4 respectively. 

The eccentricities for Au-Au collisions at yfsnjsi — 
200 GeV as function of impact parameter are shown in 
fig. [U and the corresponding plot for Pb-Pb at yfswN = 
2.76 TeV is shown in fig. [U The error bars show the 
fluctuations a £n — y/ (e 2 ) — (£ n ) 2 and, assuming not too 
large fluctu ation s in the final state evolution, the Au-Au 
ratio a F J \/ (e 2 ) agrees with data for o- V2 /v 2 from PHO- 
BOS [29| and STAR [H for all centralities. 

The angle <f> n between the participant plane of the n:th 
moment and the event plane can be seen in fig. El As seen 
in 0, 02 is strongly correlated to the event plane angle 
f i, at medium impact parameters where the geometric 
effect is strong, and significantly less correlated at central 
events where the ellipticity is more fluctuation driven. 03 
is independent of impact parameter, and cf>i only weakly 
correlated. 



IV. CORRELATIONS AND FLUCTUATIONS 

With the dynamics from the saturated BFKL cascade, 
DIPSY is expected to describe correlations and fluctua- 
tions between different rapidity slices within one event. 

One such observable is the ellipticity in the forward 
region as function of the ellipticity in the backward re- 
gion, seen in figure 0J The ellipticity is seen to be more 
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FIG. 3: (Colour online) The angle between the event plane 
and the participant plane for N par t > 100 Au-Au collisions 
at y^ijvlv" = 200 GeV. The normalisation is such that uncor- 
rected angles yield a constant 1. 



FIG. 5: (Colour online) The angle between (f> n in the forward 
region and backward region at N par t > 100 Au-Au collisions 
for y^ijvjv = 200 GeV. The normalisation is such that uncor- 
rected angles yield a constant 1. 
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FIG. 4: (Colour online) The correlation between e n in the 
forward region and backward region at N par t > 100 Au-Au 
collisions for ^/snn = 200 GeV. n = 2 (left) and n = 3 (right). 

correlated over rapidity than the triangularity. The cor- 
relation between e n in r/ E (1, 3) and in 77 G (—3, —1) can 
be quantified with the correlation coefficient 



"' £ " V((e F n) 2 )-(e F n)W((^) 2 )-(e^) 2 

where the index F and B implies the eccentricity in the 
forward, 77 £ (1,3), or backward, r/ G (—3,-1), region. 
p n is 1 if the eccentricities are perfectly correlated, and 
if they are completely uncorrelated. At N part > 100 
RHIC Au-Au, p e F e B, or p n for short, is 0.49, 0.81, 0.57, 
0.58 for n = 1, 2, 3, 4 respectively. Notice that £2 is more 
correlated over rapidity as its origin is a systematic effect. 
The same is seen at LHC Pb-Pb, where yOi,2,3,4 are 0.62, 
0.88, 0.71, 0.69 respectively. 

Looking at the angle between the participant plane in 
the forward and backward region in fig. [SJ it is seen that 
the strong correlation in the ellipticity is coming from the 
mid-centralities, where the systemtic effect is strongest. 

Experimental verification of these correlations in the 
angle and magnitude of the corresponding anisotropic 
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FIG. 6: (Colour online) The number of correlated gluons 
ricorr, and the correlated gluons typical transverse range, in 
pseudorapidity (2, 3) due to a trigger gluon in pseudorapidity 
(-3, -2) in b = Au-Au with ^/s^iv = 200 GeV. R is the 
distance from the center of the collision. 



flow observables would reinforce the picture of flow com- 
ing from systematic effects and fluctuations in the initial 
transverse geometry. 

A more direct measure of the positional correlation 
between rapidities is to measure the transverse density of 
gluons in the forward rapidity slice triggered by a gluon 
in the backward rapidity slice. An extra n corr gluons 
are seen on top of the background in the corresponding 
transverse position in the forward slice, due to the gluon 
chain from the trigger gluon that may pass through also 
the opposite rapidity bin. This correlation is analogous to 
a flux tube in other models. Fig. [5] shows central (b = 0) 
Au-Au collisions with trigger particles at a distance R 
from the center of the collision. 

It is seen that in the center of the collision, the most 
dense region, fewer gluons are correlated to the trigger 
gluon. This is a sign of saturation, as the trigger gluon 
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chain does not propagate unhindered to the other ra- 
pidity slice, but has a large probability of merging with 
other chains before that, giving fewer correlated gluons in 
the other slice. In the peripheral region of the collisions, 
n corr « 0.4, which is similar to pp collisions in DIPSY. 
This shows that the approximation of binary collisions, 
where flux tubes propagate independently through ra- 
pidity, breaks down in the dense region of a heavy ion 
collisions. 

The shorter range of the correlation in the center of 
the collision is again due to the saturated environment 
preferring smaller dipoles through the dipole swing. 

V. CONCLUSIONS 

A new method to generate t = states of high en- 
ergy heavy ion collisions has been introduced with DIPSY. 
The model has been tuned to pp and j*p minimum bias 
events, and no new parameters have been introduced for 
heavy ions. DIPSY is based on BFKL and includes all 
mergings and splittings of gluon chains, describing all 
fluctuations and correlations in a saturated environment. 

t = events for Au-Au and Cu-Au at RHIC, and Pb- 
Pb at LHC are generated and available online, see sec.llTl 

Eccentricities £1,2,3,4 and their angles ^1,2,3,4 have 
been studied and give results similar to other mod- 
els, with fluctuation-driven quantities generally slightly 
larger. Predictions were made for correlations between 
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